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ABSTRACT 

Perturbations  of  a plane  premixed  flame  propagating  steadily  along  a 
uniform  duct  (cf.  MRC  TSR  #1818)  are  considered.  These  may  be  caused  by 
heat  loss  through  the  sides  of  the  duct,  for  which  conditions  are  steady, 
or  by  slight  variations  in  the  cross  section,  when  a slowly  varying  flame 
results.  Such  unsteadiness  may  also  be  self-induced . The  results,  which 
concern  quenching,  stability  and  velocity  changes  as  well  as  the  behavior 
of  an  elementary  burner,  are  remarkably  sensitive  to  whether  the  Lewis 
number  is  greater  or  less  than  one. 
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MATHEMATICAL  THEORY  OF  IAMINAR  COMBUSTION.  III.  PERTURBATIONS 
J.  D.  Buckmaster  and  G.  S.  S.  Ludford 

1.  Introduction 

* 

The  most  striking  feature  of  the  result  (11.22)  is  that,  for  fixed  D and  Js> 
changes  in  Tm  produce  much  larger  changes  in  M because  of  the  0 in  the  exponential. 
More  precisely,  an  0(0  S change  in  flame  temperature,  due  to  such  a change  in  T 
for  the  unbounded  flame,  produces  an  0(1)  change  in  the  burning  rate.  Such  a perturba- 
tion in  flame  temperature  may  also  be  caused  by  heat  loss  through  the  sides  of  a uniform 
duct  along  which  a flame  is  propagating,  so  that  we  may  expect  a similar  change  in  its 
speed . 

Perturbations  may  also  be  produced  by  slight  variations  in  cross  section  of  the  duct, 
giving  rise  to  the  new  feature  of  unsteadiness.  Now  there  will  be  slow  variations  in 
the  combustion  field,  developing  on  a time  scale  0(0).  Indeed  such  slow  variations  can 
be  self-induced  by  residual  perturbations  of  the  initial  conditions  (on  that  time  scale) 
in  the  absence  of  boundary  perturbations.  In  all  such  cases  an  obvious  conjecture  is 
that  the  flame  velocity  is  not  close  to  the  unperturbed  value. 

It  is  not  so  obvious  that  slowly  varying  flames  behave  quite  differently  according 
as  the  Lewis  number  is  smaller  or  greater  than  one.  Moreover,  l = 1 is  exceptional ; 
in  particular  there  are  no  self-induced  variations.  In  order  to  develop  the  theory  fully, 
it  is  necessary  to  re-do  Sec.  II. 4 for  general  Lewis  numi  r,  and  that  is  where  we  start. 
The  remainder  of  the  chapter  is  an  introduction  to  slowly  varying  flames,  unearthed  by 
a perturbation  theory  which  is  proving  increasingly  useful  in  understanding  the  true  role 
of  such  idealized  solutions  as  the  premixed  plane  flame. 

2.  Modifications  for  £ * 1 

* 

The  analysis  follows  that  in  Sec.  II. 4 closely.  The  flame  sheet  is  again  located 
at  x » beyond  it  the  temperature  is  still  constant  and  there  is  no.-ieactant. 

* 

This  notation  refers  to  the  equation  or  section  number  in  Part  II,  MRC  TSR  #1818. 
Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-75-C-0024. 
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because  there  is  no  Shvab-Zeldovich  variable,  the  latter  fact  must  now  be  deduced  from 


(1)  Q(Y,f)  • - ft(T , 1 ) everywhere  . 

For  the  same  reason  both  Y a vi  T must  be  carried,  aixl  we  write 

(2)  Y - Y ♦ X~V  </X-l),  T = T + T'(e*-1)  for  0 < X < *.  , 

s S s s * 

where  Ys*Ys*rs,Tg  dre  •'‘'I1108  at  x ■ 0 determined  by  how  the  reactant  is  supplied 
there.  (Fig.  II. 1 stil!  wives  the  qeneral  shape  of  these  profiles.)  Continuity  at  the 
flame  sheet  now  shows  that 

(3)  x - X " in ( 1-X Y /Y')  - Kn 1 1 ♦ IT  -T  )/T'l  . 

■ S s s s 

The  discussion  of  reactant  supplies  in  Sec.  II. 3 showed  that  only  three  conditions 

are  provided  at  x = 0,  so  that  we  may  expect  Y’  and  T (say)  to  he  related  to 

s •*> 

Ya,Ts#T^.  The  result  (3)  provides  one  such  relation  and  the  second,  namely 

(4)  T - Y + T - (X"lY'*T‘)  » J + E 

*»  s s s s s s 

comes  from  a heat  balance  between  the  stations  x » 0 and  * [formally  by  integrating 
equation  ( 1 )).  The  second  enables  the  first  to  bo  re-written 


(5)  r’1Y'(l-^Y  /Y')1’1^  - -T*  , 

8 S S S 

so  that  now  the  reductions  (II. 6)  and  Y‘  ■ -T'  for  X = 1 are  clear. 

s s 

To  ensure  0 < xt  < ” we  must  have  Y^  < 0 and  hence  T^  > 0:  whatever  supplies 

the  reactant  must  still  be  a (conductive)  heat  sink.  Since  the  T-profile  is  monotonic, 

it  follows  that  the  temperature  of  the  flame  is  always  higher  than  that  of  the  supply, 

something  that  is  not  obvious  from  the  relations  (4)  and  (5)  when  X * 1. 

Ys  * Tg  is  the  adiabatic  flame  temperature,  so-called  because  it  is  the  final 

temperature  when  the  reactant  is  allowed  to  decompose  undor  adiabatic  conditions.  However, 

the  system  receives  heat  -X  *Y'  by  diffusion  and  loses  heat  T'  by  conduction,  where 

s s 

the  relation  (5)  shows 

(6)  -X  *Y'  - T*  according  as  X — 1 . 

s < s > 

Thus  the  actual  flame  temperature  (4)  is  higher  or  lower  than  the  adiabatic. 


. 
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Here  the  expansions  (11.17)  lead  to  the  structure  equations 

(7)  £_1d2y/d£2  - -d2t/dt2  = Aye*"  , 
which  matching  requires  to  bo  solved  under  the  boundary  conditions 

(8)  y « * o (1) , t = Jg£  ♦ o(l)  as  £ -♦  — y,t  = o(l)  as  £ +°°  » 

where 

(9)  J = Y - jf'V  . 

S S s 

now.  We  immediately  conclude  that  there  is  a local  Shvab-Zeldovich  variable  and  that  it 
otays  constant  through  the  layer,  more  precisely 

(10)  t + £-1y  = 0 . 

The  first  integral  (11.20)  is  therefore  replaced  by 

(ID  (dt/d£)  2 = 2JCA11  + ( t— 1 ) et] 

and  the  result  (IX. 22)  by 

(12)  A - DM*2  = J262exp(0/T  )/2£t4  , 

S oo  oo 

but  the  structure  (11.23)  still  holds.  The  formula  (12),  as  well  as  its  correction  of 
relative  order  0 1,  is  due  to  Bush  & Fendell  (1970). 

The  modifications  to  Secs.  1 1. 5 and  6 which  have  to  be  made  when  £ =*=  1 are  quite 
similar.  Of  more  importance  for  our  purposes  is  the  form  which  the  results  take  for  the 
unbounded  flame.  By  transferring  the  origin  to  the  flame  sheet  and  then  letting  Tg  -*•  0 
we  find 

(13)  Y = Y_oo ( 1-e  x),  T = T_m  + for  x < 0 , 

(14)  Y = 0,  T = Tb  = Y_oo  + for  x > 0 

and  the  eigenvalue  (12)  with  Js  = Y „ again.  Note  that  the  adiabatic  flame  temperature 
obtains  here  because  there  are  no  fluxes  into  or  out  of  the  system. 

Suppose  now  that  is  increased  by  fit  . Then  the  flame  temperature  changes  to 

(15)  T = T + 6t 

♦ 00  oo 

and  the  formula  (12)  shows  that 
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(lb)  M = M„e  “ where  M„  » /2jCd  T exp(-8/2T  )/J  8 

is  the  burning  rate  for  the  original  flame  temperature  T . This  result  has  universal 
validity,  i.e.  it  holds  whatever  the  nature  of  the  perturbation,  steady  or  not.  The 
reason  is  that  in  determining  the  eigenvalue,  the  perturbation  only  intrudes  through 
the  matching  of  t at  £ = +«*  (which  leads  to  the  exponential  factor)  . 

3.  Theory  of  Perturbations 

The  equations  that  we  shall  discuss  are 

(17)  3p/3t  + 3(pv)/3x  = Sf  , 

(18)  p3Y/3t  + pv3Y/3x  - •£  232Y/3x2  = -AY  exp(-0/T)  + Sg  , 

(19)  p3T/3t  + pv 3T/3x  - 32T/3x2  = AY  exp(-8/T)  + 6h  . 

Here  f,g,h  should  be  regarded,  for  the  moment,  as  arbitrary  functionals;  note  that  the 
Lewis  number  has  been  reinstated.  For  propagation  into  a quiescent  gas, 

(20)  Y , T, v -♦  Y a>,  T_qo,  0 as  x -«>  . 

The  units  by  which  these  equations  have  been  made  dimensionless  are  based  on  a representa- 
tive mass  flux,  which  will  be  taken  as  that  through  the  flame  sheet  in  the  absence  of 
perturbations,  namely  the  Mq  of  equations  (16);  so  that  A = DM~2  in  the  above  equations. 

To  avoid  an  unnecessary  factor  we  have  changed  the  unit  of  density  from  p to 

c 

P_ot)Pc  ~ p /T_^  (with  corresponding  changes  in  other  units).  The  incoming  velocity  is 
then  1 and  T_^  appears  in  the  equation  of  state,  i.e. 

(21)  p = T /T  . 

— ao 

For  finite  0,  it  is  possible  to  define  a flame  speed  only  when  the  combustion  is 
steady:  when  the  structure  changes  with  time  in  all  frames  of  reference,  the  location  of 

the  flame  (and  hence  its  speed)  is  not  a precise  concept.  In  the  limit  0 -*•  <*,  however, 
it  becomes  precise  since  the  reaction  is  confined  to  a sheet.  With  x#(t)  denoting 
the  position  of  the  sheet,  we  write 

(22)  s - x# (t)  - x 


! 

i 


.. 

I 

i 

L 
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for  the*  distance  ahead  to  obtain 


i 


(23)  3p/3t  - 3 IP  (V.v)  ]/3s  * 6f  , 

-12  2 

(24)  p3Y/'tt  - p(V+v)3Y/3s  - £ 3 Y/3s  = -AY  exp(-6/T)  ♦ 6g  , 

(25)  p3T/3t  - p (V+v) 3T/3s  - 32T/3s2  = AY  expl-fl/T)  t 6h 

where  V = x#(t)  is  the  flame  speed.  The  unper turbed  flame  corresponds  to  V = 1 and  we 
are  concerned  with  flames  for  which  V stays  0(1).  We  may  expect  that  to  be  the  case 
for  f,g,h  all  0(1)  and,  when  conditions  are  unsteady,  3/3t  = 0(6)  i.e.  slowly 

varying  flames.  Then  the  result  ( 1 6a ) , which  now  reads 

t J7 

(26)  V *=  e 


because  p is  not  changed  to  0(1),  shows  that  our  task  is  to  evaluate  the  increase 
6t.^  in  flame  temperature  due  to  the  perturbations.  Before  doinq  that,  we  note  for  future 
reference  that 


(27)  Y 0 behind  the  flame  sheet  . 

2 

If  it  were  not,  the  reaction  term  in  equation  (70)  would  be  0(0  ),  since  T = T^  to 
leading  order,  and  hence  unbalanced.  Otherwise  stated,  Y relaxes  to  zero  in  a time  0(0  ^) . 

The  first  calculation  of  t^  was  for  distributed  heat  loss,  which  was  considered 
independently  by  Buckmaster  (1976)  and  by  Joulin  & Clavin  (1976)  using  straightforward 
matching  of  3 terms  in  the  reaction  zone  with  2 terms  on  either  side,  i.e.  one  more  term 
in  each  expansion  than  is  required  without  loss.  Other  perturbat ions  were  later  treated 
by  Buckmaster  (1977),  who  by  then  realized  that  tw  (for  whose  determination  the  extra 
terms  are  introduced)  can  be  calculated  directly  from  the  overall  change  in  enthalpy  of 
the  mixture  up  to  and  including  the  flame  sheet. 

The  direct  calculation  is  as  follows.  If  equations  (24)  and  (25)  are  added  to 
eliminate  the  reaction  terms  and  then  integrated  from  s * 0-  to  **,  we  obtain 


/ 

0- 


P(v+V) 

3 s 


3_(_Y+T) 

3t 


ds  , 


correct  to  0(6),  since  Y vanishes  behind  the  flame  sheet.  Integrating  by  parts  on 
the  left-hand  side  and  using  the  continuity  equation  (23)  now  yields 
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I **  r 

(T.-T.,V  ’ H * I ~ * 0 H - (y+T-T.)5f  - 6g  - 6h]ds 

'o-  0-  9 J 

in  view  of  the  equation  of  stJue  (73).  Here  T # = T(0-,t)  may  properly  be*  called  the 
flame  temperature,  while  is  the  adiabatic  flame  temperature  Y ♦ T . Thus 


(28) 


t V 


I silo  - + £ [<T”~Y)  ^ ' 0 S * (Y+T-T^)  f + g ♦ hjds  , 


where 


(29)  T = 6t 

is  a long  time  (first  introduced  by  Sivashinsky,  1974)  and  all  terms  are  to  be  evaluated 
to  leading  order  only. 

The  whole  problem  is  reduced  to  calculating  the  right-hand  side  of  equation  (28) 

in  terms  of  V for  the  particular  perturbation  of  interest.  It  then  becomes  an  equation 

for  the  flame  speed  (the  goal  of  our  analysis)  since  the  universal  formula  (26)  shows 

2 

the  left-hand  side  to  be  V In  V . 


4.  Steady  Heat  Loss 

In  practice  the  mixture  must  be  confined  laterally  and  there  is  then  the  possibility 
of  losing  heat  to  the  sides.  For  example,  heat  conduction  through  the  wall  of  the  tube 
introduced  in  Sec.  II. 2 will  result  in  a heat  loss  proportional  to  T - T , if  we 
assume  that  the  temperature  outside  the  tube  at  any  point  is  that  of  the  mixture  at 
x « -°°.  Similarly,  thermal  radiation  (as  distinct  from  chemiluminescence)  is  present, 

though  usually  negligible,  in  combustion  processes,  resulting  in  loss  proportional  to 

4 4 

T - T oo.  The  importance  of  considering  such  effects  comes  from  the  observed  inability 
of  the  mixture  to  ignite  if  the  tube  is  too  thin,  a phenomenon  about  which  our  analysis 
so  far  has  nothing  to  say.  Then  conductive  loss,  though  not  radiative,  is  significant. 

Heat  losses  are  normally  minimized  as  much  as  possible,  so  that  the  propor t ional i ty 
factors  are  small.  Our  analysis  will  therefore  be  based  on  0(9  *)  factors,  but  the 
dependence  of  the  loss  on  temperature  will  be  quite  general;  in  particular,  both  conductive 
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and  radiative  losses  are  covered.  A term  -i$M  2<t(T)  is  added  to  the  right-hand  side  of 
the  temperature  equation  (II. 4),  where  $ vanishes  at  T but  is  otherwise  an 


- 2 

arbitrary  positive  function,  and  6 has  the  definition  (11.16) . (The  factor  M ensures 


that  $ is  independent  of  the  mass-flux  and  density  units.)  In  particular. 


(30) 


k (T-T  ) for  conductive  losses 


4 4 

k(T  -T  ) for  radiative  losses 


where  k is  a given  constant. 

We  are  therefore  dealing  with  the  special  case  of  equations  (23-25)  in  which 


(31) 


-2. 


3/3t  = f = g = 0,  h = -Mq  ifi  (T)  . 


The  continuity  equation  integrates  to 

(32)  o(V+v)  = V , 

with  V constant,  so  that  to  leading  order  we  are  left  with  the  unperturbed  steady 
equations.  Use  can  therefore  be  made  of  the  result  (13)  to  write 


(33) 


—£  v s —V 

Y = Y^U-e  ),  T = T-oo  + Y^e  for  s > 0 , 


but  (14)  must  be  taken  one  term  further.  From  the  perturbed  temperature  equation  we  easily 
f ind 


(34)  T = T^  + 6 [ + V_1<MTJs]  for  s < 0 , 

a result  that  is  needed  to  calculate  3t/3s  at  0-.  [An  exponentially  growing  term  has 
been  discarded  as  unmatchable  with  the  decay  to  ambient  temperature  which  occurs  on  a 
scale  0(6)  . ] 

The  right-hand  side  of  (28)  may  now  be  calculated,  giving 


(35) 


t = -tv  2 where  = t (T  ) + / 4(T  +Y  e °)do  , 

OD  0 00  J —or.  —on 

0 


—00  —00 


V2fn  V2  + t = 0 


so  that 
(36) 

is  the  equation  for  the  flame  speed.  The  result  is  graphed  in  Fig.  ]j  such  C-shaped 
response  curves  are  conmon  in  combustion  theory.  The  function  t consists  of  two  terms: 


i 

I 
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Fiq.  1:  Flame  velocity  V versus  heat-loss  parameter 

The  notation  v^tO)  and  #(0)  refers  to  Sec.  7 

2 

where  S,  the  point  at  which  #(0)  = 4v_^(0) /3 , 
plays  a role. 
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tile  integral  represents  the  total  heat  loss  up  to  the  flame  while  the  other  accounts  for 
conduction  into  the  burnt  mixture  behind. 

The  top  branch  hows  that  heat  loss  decreases  the  flame  velocity,  though  if  f 
exceeds  the  critical  value  e ' there  is  no  corresponding  V,  i.e.  no  flame.  The  term 

extinction  is  therefore  appropriate  for  the  point  at  thich  the  curve  turns  back,  the 

-1/2 

burning  rate  then  being  e ( = 0.61)  times  the  adiabatic  value  for  the  same  pressure. 

(Note  that  this  result  is  independent  of  all  parameters.)  The  bottom  branch  is  generally 
believed  to  be  unstable  (Qnmons,  1071)  and  hence  of  no  interest,  but  we  shall  postpone 
a discussion  of  that  point  until  the  next  section. 

It  has  been  known  for  over  150  years  that  a flame  will  not  propagate  through  a wire 
gauze,  the  principle  underlying  the  miner’s  safety  lamp  invented  by  Davy.  This  is  suggested 
by  the  present  results:  a cold  wire  gauze  of  sufficiently  small  mesh  is  an  effective 
heat  sink  and  corresponds  to  values  of  $ greater  than  critical. 

Extinction  in  small  tubes  can  now  be  explained.  Heat  loss  by  conduction  through  the 
wall  is  propxirtional  to  surface  area  and  hence  the  radius,  while  the  volume  of  the 
mixture  is  proportional  to  its  square.  The  heat  loss  per  unit  volume  therefore  varies 
inversely  as  the  radius,  and  hence  so  does  the  parameter  $ . There  is  therefore  a 
definite  tube  radius,  independent  of  pressure,  below  which  steady  combustion  will  not  occur. 

5 . Slowly  Varying  Flames 

Flames  can  be  unsteady  for  a variety  of  reasons.  Evolution  from  ignition,  passage 
down  a tube  of  varying  cross  section  and  instability  are  all  examples  of  unsteady  combus- 
tion. Here  wc  are  concerned  with  the  one-dimensional  problem  of  a flame  evolving  from 
initial  data;  in  general,  such  an  evolution  can  be  expected  to  take  place  on  an  0(1) 
time  scale.  Indeed  there  have  been  numerical  integrations  (with  f,g,h  all  zero)  that 
show  the  emergence  of  what  appears  to  be  a steady  wave  from  a mass  of  hot  quiescent 
gas  on  such  a time  scale  (Zeldovich  & Barenblatt,  1959) . 

Our  techniques  cannot  describe  such  phenomena  in  general  since  the  right  side  (28) 
will  not  be  small.  However,  from  the  initial  development  there  emerges  a structure, 
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characterized  by  the  single  variable  V,  that  undergoes  a slower  evolution  on  the  time 
scale  0,  and  during  this  phase  jf  the  unsteady  process  the  time  derivatives  are  perturba- 
tions that  can  be  incorporated  into  the  analysis  as  easily  as  heat  loss,  for  example. 

To  leading  order  the  continuity  equation  still  has  the  integral  (32),  where  V is 
now  a function  of  i,  sc  that  the  profiles  (33)  remain  valid  ahead  of  the  flam«-  sheet. 

To  them  we  must  add 

(37)  p - T yiT^iY  _>e~VS) 

but  an  expression  for  v w*ll  not  be  needed.  As  before,  the  remaining  question  is  to 
calculate  3T/9a  .it  ()-,  whicn  comes  from  the  temperature  perturbation  behind  the 


flame  sheet. 

To  begin  with,  consider  perturbations  duo  solely  to  the  unsteadiness.  Then 
f * g - h • 0 and  hence 

(3H)  T ■ To  + fit  for  s < 0 , 

so  that 

1 jC-1 


( 3‘>) 


t -bV  *V  where  b ■ Y T / ° — do 

* -a,  -a,  J *r  * Y Cl 

U — u>  —a) 


and  the  dot  denotes  rate  of  change  with  respect  to  the  time  i.  The  equation  for  V(i) 
is  therefore 


(40) 


bV  + V3?n  V"’  - 0 . 


Note  that  the  parameter  b changes  from  jositivu  to  negative  as  the  Lewis  number  increases 
through  one,  cjiving  the  first  indication  that  f - 1 is  not  typical. 

We  are  only  free  to  specify  the  initial  value  of  V and  then  the  present  approx ima- 
t ion  to  the  combustion  field  is  uniquely  determined.  As  mentioned  previously,  the  present 
approach  does  not  solve  the  general  ini tial -value  problem  since  the  initial  evolution 
takes  place  on  the  0(1)  time  scale.  But  f rom  this  initial  development  emerges  a 
structure,  characterized  by  the  simile  variable  V,  that  changes  on  the  0(0)  time  scale. 

When  £ > 1,  b is  negative,  anti  when  £ < l,  b is  jositive.  The  closer  £ is  to  1 
the  more  rapidly  V changes,  so  that  when  f - 1 » 0(1/0)  or  smaller  there  are  no  0(1) 
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effects  on  the  slow-time  scale  and  the  result  (36)  implies 
(41)  V = 0 or  1 . 

We  infer  that  the  initial  evolution  ends  at  one  of  these  values,  though  wc*  may  expect 
the  ultimate  value  to  be  reached  more  slowly  in  t as  the  Lewis  number  gets  further 
from  1 on  the  scale  of  1/0.  It  is  tempting  to  take  unit  Lewis  number  in  combustion 
problems  because  of  the  simplifications  stemming  from  qlobal  Shvab-Ze ldovich  variables 
(Sec.  1.7);  but  in  view  of  the  present  atypical  behavior  such  temptations  should  certainly 
be  resisted  for  unsteady  problems. 

The  values  (41)  are  also  possible  steady  states  as  t -*■  However,  the  first 
corresponds  to  t^  = <«,  when  our  analysis  breaks  down  but  indicates  that  the  flame 
temperature  differs  from  its  adiabatic  value  by  an  0(1)  amount.  A discussion  of  the 
non-uniformity,  in  the  context  of  vanishingly  small  heat  loss,  has  been  given  by 
Buckmaster  (1976).  The  second  is  unstable  when  b is  negative,  i.e.  £ > 1.  For, 
according  to  (40),  when  V is  less  than  1 it  will  be  driven  towards  0 and  when  V 
is  greater  than  1 it  will  be  driven  towards  <». 

It  might  be  thought,  therefore,  that  steady  flames  of  Lewis  number  greater  than  1 
ire  of  no  practical  interest;  but  laboratory  experience  shows  otherwise.  Moreover , if 
flames  with  f > 1 are  to  be  discarded  on  these  grounds  then  so  must  flames  with  £ \, 

since  they  arc  also  unstable  when  three-dimensional  disturbances  are  permitted  (Ch.  VI). 

The  fact  that  flames  are  stabilized  by  burners  suggests  that  velocity  and  thermal  gradients 
play  an  imi>ortant  role.  Indeed,  it  is  not  difficult  to  devise  a one-dimensional  model 
(Sec.  10)  that  stabilizes  flames  in  this  way.  However,  it  is  unlikely  that  such  mechanisms 
will  explain  stable  flames  propagating  down  uniform  tubes. 

To  be  sure,  Lewis  numbers  are  usually  not  far  from  1 so  that,  in  view  of  the 
irrelevance  of  the  slow-time  scale  when  £ - 1 =*  0(1/0),  it  is  conceivable  that  predic- 
tions of  a slow-time  analysis  are  often  not  applicable.  Alternatively,  stable  flames 
may  be  a phenomenon  of  finite  activation  energy  about  which  asymptotic  theory  provides 
faulty  information.  Finally,  stable  flames  may  be  the  result  of  chemical  kinetics  more 
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complicated  than  that  adopted  here.  Whatever  the  explanation,  activation-energy 

asymptotics  for  •£  * 1 can  provide  useful  qualitative  insights  even  though  the  appl icabi 1 i ty 

of  the  model  cannot  be  justified  at  the  present  time. 

To  determine  the  effect  of  heat  losses  we  still  let  f = g = 0 but  now  take 
-2 

h = -M0  (T)  . The  effects  of  unsteadiness  and  heat  loss  are  then  seen  to  be  additive,  i.e. 

(42)  t - -bV_2V  - »v”2  , 

00 

so  that  we  end  with 

(43)  bV  + v3en  V2  + 4>V  = 0 

as  the  equation  governing  V.  To  the  possible  steady  state  (87),  plotted  in  Fig.  1, 
roust  now  be  added 

V =»  0 for  all  $ . 

However,  the  latter  corresponds  to  t = “,  a non-uniformity  which  has  yet  to  be  treated. 

The  question  of  which  branch  of  the  multivalued  response  in  Fiq.  1 will  be  observed 
in  practice  can  be  answered  in  the  context  of  our  present  model.  The  arrows  in  Fig.  1 
give  the  direction  in  which  V changes  according  to  the  differential  equation  (43)  when 
b is  positive  (i.e.  £ < 1).  The  lower  branch  of  (36)  is  therefore  predicted  to  be 

unstable,  a result  in  accord  with  general  belief.  However,  Spalding  & Yumlu  (1959)  claim 
to  have  observed  the  slow  branch  which,  in  the  absence  of  other  effects,  can  only  be 
reconcile)  with  our  analysis  if  b is  negative  (i.e.  £ > 1)  when  the  arrows  are 

reversed.  [The  fact  that  they  used  a special  stabilizing  apparatus  suggests  that  there 
were  other  effects,  cf.  Sec.  7.J 

6.  Non-Uniform  Ducts 

We  consider  ducts  whose  cross-sectional  area  A varies  over  distances  of  order  9,  i.e. 

T 

(44)  A = <7(6x)  * <7(/  Vdi-6S)  . 

0 

In  the  usual  one-dimensional  treatment  of  such  slowly  varying  ducts,  the  divergences  in 
the  basic  equations  (1.53,54,56)  now  give  rise  to  extra  terms 
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which  correspond  to 


1 

3A 

-X'1  I 

3A 

1 ~ 

1 3A 

3T 

A 

7T  PV  , 

3s 

, A 

3s 

3s 

' A 3s 

3s 

(45)  f = F(X)pv,  g = f_1F(X)3Y/3s,  h = F(X)3T/3s  , 
where 

T 

(46)  F (X)  = O'  (X)/67(X)  with  X = / VdT  . 

0 

We  have  neglected  6s  in  the  argument  of  O for  s = 0(1),  i.e.  F is  the  logarithmic 
derivative  of  the  area  at  the  location  X of  the  flame  sheet. 

It  is  easily  seen  that  3T/3s  =0  at  0-,  so  that  the  extra  terms 

<30 

(47)  F (X)  / [£-13Y/3s  + 3T/3s  + (Y+T-TJ  (l-p)V]ds  = -bF(X) 

0 

in  t^V  come  from  the  integral . They  lead  to 

(48)  bV  + v3tn  V2-bF(X)V2  = 0 
as  the  replacement  for  the  governing  equation  (91) . 

An  immediate  consequence  is  that  a flame  sheet  propagating  steadily  (V  = 1)  down 
a duct  of  constant  cross  section  will  in  general  accelerate  on  entering  a diverging 
section  (O'  > 0)  whereas  if  it  enters  a converging  section  ( O'  < 0)  it  will  slow 
down  (Sivashinsky , 1974).  The  only  exception  is  b = 0 (X  =1),  for  then  V = 1 irrespec- 
tive of  area  changes.  Again,  f = 1 is  exceptional. 

The  F (X)  makes  the  differential  equation  (48)  very  difficult  to  discuss  in  general. 
However,  if  the  area  changes  very  slowly  indeed,  more  precisely  O(x)  = ^(ex)  where 
0 < e <<  1,  there  is  an  approximate  solution 

(49)  V = 1 + ebF(eX),  F = O' /O 

in  which  V differs  just  a little  from  1.  It  represents  a balance  between  the  effect 
of  area  change  and  that  of  enthalpy  deficiency  at  the  flame,  temporal  variations  playing 
no  role. 

There  have  been  attempts  to  discuss  the  effect  of  curvature  on  unsteady  flames  by 
supposing  that  the  difference  between  the  speed  of  a slightly  curved  flame  and  a plane 
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flame  is  proportional  to  the  cruvature,  with  the  constant  of  propor tional i ty  depending 
only  on  the  mixture  properties  ar  1 not  on  the  particular  unsteady  process.  Markstein 
(1964,  p.  22)  has  discussed  fl  .me  stability  in  this  way.  The  result  (49)  provides  a 
rational  justification  for  • uch  an  assumption,  when  it  is  realized  that  curvature  manifests 
itself  through  fractional  changes  in  area.  In  particular,  it  shows  that  the  constant  of 
proportionality  changes  sign  as  jC  passes  through  1. 

7.  An  El enentary  Flame  Ho lder 

An  extension  of  the  present  discussion  provides  some  insight  into  the  nature  of  flame 
holders  or  burners.  The  basic  function  of  a burner  is  to  maintain  a flame  in  a stable 
rest  position:  the  flame  is  brought  to  rest  by  applying  a counterflow;  stability  is 
ensured  by  appropriate  velocity  and  thermal  gradients.  Velocity  gradients  can  be  produced 
by  expansion  of  the  mixture  as  it  leaves  a tube  or  by  flow  to  the  rear  of  a bluff  body, 
while  thermal  gradients  arise  quite  naturally  by  conduction  to  the  burner,  to  mention 
some  common  examples.  These  can  be  modelled  by  a one-dimensional  formulation  of  propaga- 
tion along  a non-uniform  duct  with  a counterflow,  incorporating  heat  losses  which  depend 
on  location  as  well  as  temperature.  For  simplicity  we  shall  only  consider  conductive 
losses  through  the  side  walls  for  which  the  linear  law  (30a)  holds  with  k a slowly 
varying  function  of  x,  just  as  A is.  The  perturbation  terms  are  now 

(50)  • f = F (X)  Pv , g = £-1F(X)  3Y/3s,  h = F(X)3T/3s  - M"2k(X)(T-T  ) . 

0 — 00 

The  counterflow  changes  the  velocity  condition  ahead  of  the  flame  to 

(51)  v •*  v as  x -*•  -°° 
so  that  continuity  now  requires 

(52)  P (v+V)  - v_m  + V = 1r  (say) 

to  leading  order.  At  various  places  in  the  preceding  analysis  we  must  now  replace  V 
by  If.  As  a consequence  the  expansions  (27),  (33),  (34),  (37)  and  the  enthalpy  result 
(42)  with  the  extra  term  (47)  still  apply  if  V is  everywhere  replaced  by  If  except 
in  (l-p)V,  which  must  be  replaced  by  V - pv.  Noting  that 
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(53) 


♦ = k (X)  (T-T _J 


then  gives 

(54)  t v.  -b‘<~21r  - 4(X)V_1  + bF(X)V^-1  where  M^4(X)  = 2Y  k(X)  , 
so  that  the  governing  equation  becomes 

(55)  b V + “i^in  IT2  + 4 (X)V  - bF(X)V7/  = 0 . 

Without  counterflow,  i.e.  v = 0 and  V = V,  we  reach  a combination  of  the  two 
equations  (43)  and  (48)  as  expected. 

In  the  absence  of  heat  losses  ( 4> = 0 ) there  is  a stationary  solution 
V = 0 provided  V = 1 i.e.  v = 1 . 

The  flame  speed  relative  to  the  fresh  mixture  is  the  same  as  that  for  a uniform  duct. 

That  is,  the  area  changes  do  not  affect  the  flame  speed  to  leading  order,  in  striking 
contrast  to  the  unsteady  problem  discussed  earlier  where  there  was  no  counterflow.  Of 
course  the  difference  is  that  the  flame  sheet  itself  does  not  encounter  area  changes  here. 
The  corresponding  condition  is 

(56)  vL»®'n  v2«o  + *(0)  =0 

when  heat  losses  are  present.  Not  surprisingly,  this  is  just  the  relation  (36)  with  V 
replaced  by  v_^  and  with  $ given  its  value  at  the  origin  (i.e.  any  finite  x) ; so 
that  Fig.  1 applies.  In  either  case,  with  or  without  heat  loss,  the  location  of  the  flame 
is  indeterminate. 

When  the  flame  is  moving,  v m becomes  a function  of  X.  For  it  is  the  mass  flux 
of  the  counterflow  which  is  held  fixed,  whereas  the  area  <7(x)  at  the  flame  varies  with 
its  position.  That  is  the  area  which  is  effectively  constant  in  the  integration  (52) 
of  the  equation  of  continuity  on  the  x-scale,  so  that 

(57)  v^  « <7-1(x)  . 

Accordingly,  we  should  write  v (0)  in  the  condition  (56). 

To  investigate  the  stability  of  a stationary  flame,  time  dependent  solutions  of 
equation  (55)  are  sought  in  which  V is  small  (though  still  large  compared  to  6),  so 
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that  If  differs  just  a little  from  v The  linear  equation  governing  small  changes 

from  the  rest  position  is 

(58)  X + (b-1D  - 2F(0)  Jv^fOJX  + |$'(0)  - DF(0)v  ^ (0)  ] b"’  v_^  (0)  X = 0 , 
where 

(59)  D = 2v  (0)11  + 2 in  v (0) ) = -d$(0)/dv  (0)  ; 

— oo  — OO  —Ob 

and  the  necessary  and  sufficient  condition  for  stability  is  that  the  coefficients  of  both 
X and  X are  positive.  The  stability  criteria  are  therefore 

(60)  2F(0)  < b and  4>'(0)  J DF(0)v  ^(0)  according  as  b J 0 . 

To  see  how  these  conditions  can  be  met,  consider  a duct  whose  cross  section  varies 
in  size  but  not  in  shape.  The  argument  used  for  a circular  tube  at  the  end  of  Sec.  4 
then  shows  that  the  heat-loss  coefficient  varies  inversely  as  a linear  dimension  of  the 
cross  section,  i.e. 

(61)  k =<7'1/2(X)  . 

It  follows  that  t'(0)  = - ^ 4 (0)  F (0)  , so  that  the  second  criterion  (60)  may  be  written 

(62)  ♦ (0) F (0)  £ j v^(0)F(0)  according  as  b J 0 , 
since  Dv^lO)  = 2(^(0)  - *(0)). 

To  interpret  these  conditions,  note  that  the  derivative  D is  positive/negative  on 
the  upper/lower  branches  in  Fig.  1.  When  the  Lewis  number  is  less  than  1,  i.e.  b > 0, 
a flame  is  stable  only  if  it  corresponds  to  a point  on  the  curve  above  S and  the  duct 
converges;  indeed  on  the  lower  branch  the  convergence  must  be  sufficiently  rapid,  i.e. 

(63)  F (0)  < D/2b  . 

When  the  Lewis  number  exceeds  1,  i.e.  b < 0,  there  is  stability  only  on  the  lower 
branch  and  then  the  duct  must  converge  or,  for  points  above  S,  diverge  sufficiently  slowly: 

(64)  F (0)  < D/2b  . 

A well-insulated  duct  must  therefore  converge  and  then  the  flame  with  v_^  small 
is  stable  for  X > l while  that  with  v_a  close  to  1 is  stable  for  X < l.  The  results 
for  a uniform  duct  in  Sec.  5 are  a limiting  form  of  these  conclusions. 
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